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Abstract 

The results of a numerical study to assess the effect 
of skin buckling on the internal load distribution in a 
stiffened fuselage panel, with and without longitudinal 
cracks, are presented. In addition, the impact of changes 
in the internal loads on the fatigue life and residual 
strength of a fuselage panel is assessed. A generic nar- 
row-body fuselage panel is considered. The entire panel 
is modeled using shell elements and considerable detail 
is included to represent the geometric-nonlinear re- 
sponse of the buckled skin, cross section deformation of 
the stiffening components, and details of the skin-string- 
er attachment with discrete fasteners. Results are pre- 
sented for a fixed internal pressure and various 
combinations of axial tension or compression loads. Re- 
sults illustrating the effect of skin buckling on the stress 
distribution in the skin and stringer, and fastener loads 
are presented. Results are presented for the pristine struc- 
ture, and for cases where damage is introduced in the 
form of a longitudinal crack adjacent to the stringer, or 
failed fastener elements. The results indicate that axial 
compression loads and skin buckling can have a signifi- 
cant effect on the circumferential stress in the skin, and 
fastener loads, which will influence damage initiation, 
and a comparable effect on stress intensity factors for 
cases with cracks. The effects on stress intensity factors 
will influence damage propagation rates and the residual 
strength of the panel. 

Introduction 

Buckling of a fuselage panel causes local displace- 
ment and stress gradients that are greater than the gradi- 
ents in a panel that is not buckled. These local gradients 
influence the failure-initiation, damage-propagation and 
residual-strength behavior of a fuselage shell. Thus, if a 
fuselage shell is allowed to buckle, then these gradients 
must be understood in order to predict accurately fuse- 
lage structural integrity and residual strength. 
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Transport-fuselage shells are designed to support 
internal pressure and mechanical flight loads which result 
in local panel loads that consist of various levels of longi- 
tudinal tension or compression, circumferential tension, 
and shear, and that cause a geometrically nonlinear shell 
response. Typical metallic fuselage structure consists of 
built-up stiffened panels with a thin skin attached Co lon- 
gitudinal stringers and circumferential frames. To maxi- 
mize structural efficiency, fuselage shells are usually 
designed to allow the fuselage skin to buckle in compres- 
sion and shear at a load level that is below the design limit 
load condition for the shell. Thus, it is assumed that 
cracks may exist in the structure during the service life of 
the aircraft, and that loading conditions could occur that 
would cause the fuselage skin with cracks to buckle. 

Skin buckling causes nonlinear deformations and 
changes in the stress distribution in the skin, the internal 
structure, and the fasteners connecting the skin and the in- 
ternal structure. Failure initiation and propagation in the 
built-up structure may involve crack initiation in the skin 
or stiffening elements, or fatigue or strength failure of the 
fastener elements connecting the skin to the stiffening el- 
ements. The structural integrity of a built-up structure 
subjected to combinations of internal pressure and me- 
chanical Right loads can be studied analytically with a 
nonlinear structural analysis capability, but a high-fideli- 
ty modeling and analysis methodology must be applied to 
obtain accurate predictions of the state of stress in each 
component of the structure. Most residual-strength anal- 
ysis studies reported in the literature for fuselage shells 
with cracks 1-9 have been limited to internal pressure 
loads only, where the shell is in biaxial tension. The re- 
sults of analytical studies of the nonlinear response of un- 
stiffened aluminum shells with longitudinal cracks and 
subjected to internal pressure and axial compression 
loads 10-12 have indicated that the crack-growth character- 
istics of longitudinal cracks are influenced by the biaxi- 
al-loading ratio, x, defined as the ratio of the longitudinal 
stress resultant to the circumferential stress resultant. The 
influence of biaxial loading on cracked stiffened panels 
was reported in Refs. 13 and 14, hut skin buckling was 
not considered in either reference. In addition, results of 
a fatigue test of an A300B fuselage 1 ^ indicated that com- 
pressive stress directioned parallel to a crack may in- 
crease the stress intensity factor by 40 %. 
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The present paper describes the results of a numer- 
ical study of the nonlinear response of a typical aircraft 
fuselage panel subjected to combinations of internal 
pressure and axial compression loads, and with and with- 
out a longitudinal crack. The paper focuses on the re- 
sponse of the skin and stringer region between the frames 
and describes the modeling and analysis methodology 
required to represent the nonlinear postbuckling re- 
sponse of the built-up structure with rivets attaching the 
skin to the stringers. Numerical results are presented 
first for a pristine structure, to demonstrate the effect of 
skin buckling on the fastener loads and the stress distri- 
bution in the fuselage panel, and to assess the impact of 
skin buckling on damage initiation. Then, results are 
presented for fuselage panels with cracks to assess the 
impact of skin bucking on the tendency for damage prop- 
agation. The numerical results are intended to provide 
insight into the complex nonlinear response of a built-up 
structure and to demonstrate the types ol structural re- 
sponse that can be simulated by using modern analysis 
tools. Results are presented for a generic narrow-body 
fuselage panel, and it is recognized that the quantitative 
results are sensitive to the specific configuration consid- 
ered. To allow for a qualitative measure of the effect of 
skin buckling on the response, numerical results are nor- 
malized by the response of the non-buckled structure 
subjected to nominal pressure loads. It is also recognized 
that prediction of damage initiation and damage propa- 
gation in a built-up structure typically requires a glo- 
bal-local approach, where a global shell model may be 
used to determine the structural-level response and de- 
fine the appropriate gross loads to be applied in a refined 
local model. 16 The local model may require three-di- 
mensional analysis and should be capable of represent- 
ing the effects of fastener countersink details, 
interference fit, and through-thickness crack profiles. 
The modeling and analysis effort for the present paper 
consists of a global shell model only, but considerable 
detail was included to represent the geometric nonlinear- 
ity and spatial variation in the response, as well as the 
structural-level interaction between the components in 
the built-up structure. The results from the analyses con- 
ducted for the present paper could be used as input to 
more-refined local analyses. One product from the 
present paper is an appreciation for the complexity of the 
response demonstrated by the ‘global’ model, and a bet- 
ter understanding of the types of loading conditions that 
a local model may need to represent. For the purpose of 
the present paper, estimates of the effect of skin buckling 
on damage initiation and propagation are based on stress 
results from the global shell model, and stress intensity 
factors for cracks modeled as simple, straight, 
through-cracks. 


Fuselage Panel Geometry and Finite Element Model 

The structural configuration considered is shown in 
Fig. 1, and is a generic narrow-body fuselage panel. It is 
constructed entirely of 2024-T3 aluminum alloy, with a 
74.-in. skin radius, a 0.040-in. skin thickness, Z-stringers 
with an 8. -in spacing, and Z-frames with a 20. -in. spac- 
ing. A finite element model of the stiffened fuselage 
panel with two frame-to- frame longitudinal skin bays 
and five circumferential stringer bays is shown in Fig. 1 . 
The origin of the (x,y) coordinate system shown in Fig. I 
is located on the center stringer, and midway between the 
frames. The model was defined to include one half of a 
skin-bay beyond the last stiffening member on each edge 
of the panel. The circumferential edges of the skin and 
frames have symmetry boundary conditions. The longi- 
tudinal edges of the panel have the rotational constraints 
of a line of symmetry and multi-point constraints to en- 
force a uniform longitudinal edge displacement. The 
Young’s modulus, E> for the aluminum alloy is equal to 
10.5 msi and Poisson’s ratio, v, is equal to 0.33. 

The loading condition for the fuselage panel con- 
sists of an applied internal pressure, p (which generates a 
circumferential stress resultant reaction, 7V v ), and an axial 
stress resultant, N x , which is the sum of the bulkhead 
pressure load, and an applied mechanical load. The 
stress resultants, N x and N v , represent the average load in 
pounds per inch along the longitudinal and circumferen 
tial edges of the panel, respectively. A biaxial loading 
ratio, % * is defined as the ratio of the axial load to the cir- 
cumferential load, X " N. t /Wy . A biaxial loading ratio 
X~ 0.5 corresponds to the internal-pressure-only load- 
ing condition. 

The stringer and frame elements of the panel are 
constructed of 0.050-in. -thick 2024-T3 aluminum alloy. 
The Z-frames are connected to the skin by L-shaped clips 
that have cutouts for the stringers to pass through, and a 
continuous circumferential tear strap is located between 
the L-shaped clips and the skin. The cross section di- 
mensions for the Z-stringer, Z-frame, and L-clip are 
shown in Fig. 1. In the finite-element model, the cross 
sections are modeled as branched shells, rather than dis- 
crete beams, so that cross-section deformation could be 
represented. The portions of the skin-Z-frame-L-clip 
cross section that have multiple layers are represented by 
a single plate element with multiple layers, and thus the 
layers are attached continuously, as if the layers were 
bonded. Conversely, the skin-stringer connection is 
modeled as a riveted connection by stacking finite ele- 
ment plate elements and connecting the plate elements 
with springs in a manner similar to that described in Ref. 
17. The skin-stringer connection was defined to repre- 
sent a single row of 3/1 6-in. -diameter rivets with a 1 .-in. 
pitch. The details of the skin-to-stringer attachment 
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Axial load: N x = ^N y (% = - 1 .2 to +2) 


1 .0 in. 



I— Tear strap, t = 0.05 in. 


Figure 1 . Geometry and finite element model of stiffened fuselage shell. 

model are shown in Fig. 2. The skin and stringer are rep- stiffened locally. The fastener element used in the anal- 

resented by offset layers of plate elements. A six-spring ysis consists of two rigid links connected by six springs, 

fastener element connects a single node in the skin to a These springs represent the stiffnesses associated with 

single node in the stringer, and beam elements are used extension, shearing, bending and twisting of the fastener, 

to distribute the fastener load over a 0.25 in.-by-0.25 in. The shear stiffness of the fastener (164e3 Ib/in.) was 

area in each layer. The beam elements used to distribute computed using the empirical relation defined in Ref. 1 8, 

the fastener load to the skin are defined so that they are evaluated for a 0. 1875-in. -diameter aluminum rivet join- 

stiff relative to the fastener element, but soft relative to ing a 2024-T3 aluminum alloy skin and stringers of 

the plate elements, so that the skin and stringer are not thicknesses 0.040 in. and 0.050 in., respectively. Theex- 
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Contact elements 

Six-spring fastener element 

Beam elements to distribute 
fastener load over a 
0.25-in. -by -0.25-in. area 



Symmetric membrane 
load, K f 



Anti-symmetric 
membrane load, K u 



Symmetric bending 
load, kj 



Anti-symmetric 
bending load, k 2 


Figure 2. Detail of skin-to-stringer attachment 
model. 


Figure 3. Stress intensity factors for symmetric 
and anti-symmetric membrane and 
bending loads. 


tensional. flexural, and torsional stiffnesses of the fasten- 
er element were computed assuming that the rivet 
behaves like a simple 0. 1 875-in. -diameter aluminum 
rod. Pad-contact elements 19 are defined between the 
skin and stringers to represent the contact that would oc- 
cur if the layers of elements were to move toward each 
other. 

The finite element discretization was defined to 
represent adequately the buckled skin response (at least 
four elements per half wave). Mesh-transition elements 
which provide 2: 1 mesh refinement were used to refine 
the mesh near the stringers, particularly at the center 
stringer, where cracks are to be introduced. The result- 
ing finite element mesh had 1 57,000 degrees of freedom. 

Analysis Procedure 


The STAGS (STructural Analysis of General 
Shells) nonlinear shell analysis code, 19 which has special 
features for modeling fastener elements, contact between 
built-up components, and cracks in shell structures, was 
used to conduct the analyses. STAGS is a finite element 
code that can include the effects of geometric and mate- 
rial nonlinearities in the analysis. The code uses both the 
modified and full Newton methods for its nonlinear solu- 
tion algorithms, and accounts for large rotations in a shell 
by using a co-rotational algorithm at the element level. 
The Riks pseudo arc-length path following method 20, 21 
is used to continue a solution past the limit points of a 
nonlinear response. In addition to the nonlinear equilib- 
rium solution, output from STAGS calculations include 
the following crack-tip fracture parameters: strain-ener- 
gy-release rates (used to predict crack growth rates for 
fatigue loading conditions and residual strength crack 
extension from an elastic analysis) and the crack-lip 
opening angle (CTOA, used to determine residual 


strength crack extension from an elastic-plastic analy- 
sis). Procedures for computing the strain-energy-release 
rates using the modified crack closure integral method, 
and computing the linear-elastic stress-intensity factors 
from the strain-energy-release rates are described in 
Refs. 22 and 23. Stress intensity factors for symmetric 
and anti-symmetric membrane (Kj, K //) and bending {k h 
k 2 ) modes, shown schematically in Fig. 3, can be com- 
puted from individual components of the strain-ener 
gy-release rate. For the present study, the total stress 
intensity factor K T is calculated from the total strain-en- 
ergy-release rate, G : 


k t = Jeg (i) 

To provide a qualitative measure of the effect of skin 
buckling on the response, stress intensity factor results 
in the current paper (K h *//- k h k 2 , Kj) are defined as 
the stress intensity factors described above, normalized 
by Kfiox a biaxial loading ratio value x - 0.5. 

The STAGS code allows for separation of the ap- 
plied loads into two load sets that can be scaled indepen- 
dently during a nonlinear analysis. To conduct the 
analyses for the current paper, an internal pressure was 
applied in one load set, while an axial load was applied 
in the second load set. One nonlinear analysis was con- 
ducted to obtain the solution for an internal pressure of 8 
psi, and a large axial tension corresponding to a biaxial 
loading ratio x= 2.0 , or four times the bulkhead tension 
load for a pressure-only loading condition. Then, a sec- 
ond nonlinear analysis was conducted, keeping the pres- 
sure at 8 psi, and reducing the axial tension load to zero 
and then applying axial compression loads. The maxi- 
mum compression load considered corresponded to a bi- 
axial loading ratio -1*2, which was two times the 
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axial compressive load required to buckle the skin of the 
fuselage panel. All analyses were conducted using a lin- 
ear-elastic material model. While post-processing the 
solutions, a yield stress of 50 ksi was assumed to assess 
the validity of the elastic solutions. 

Assessment of Damage Initiation and Propagation 

Paris’ fatigue crack growth law 24 is applied to esti- 
mate the effect of changes in stress levels and stress in- 
tensity factors on damage initiation and propagation. 
Paris’ crack growth law is 


Tn = Ci&K) 


( 2 ) 


where C and // are material dependent variables which 
are determined by experiment. For aluminum alloys, 
n = 4 is a reasonable estimate. If the loading cycles are 
from zero to maximum load, then AK = K. For small 
cracks, K is linearly proportional to the remote stress 
perpendicular to the crack, and for crack initiation, the 
initial crack size is related to micro-structural disconti- 
nuities. For a given configuration subjected to two dif- 
ferent loading conditions, Paris’s law, with n = 4, states 
that the crack growth rate for each loading condition is 
proportional to the remote stress or stress intensity fac- 
tor, raised to the fourth power. The number of cycles 
required to extend a crack a distance da for some value 
of x, normalized by the number of cycles to extend the 
crack the same distance for x = 0.5, is 

N = (K 0 ,/Af x ) 4 = (o 05 /a x ) 4 (3) 

It is assumed that for mixed-mode loading conditions, 
the crack growth rate is a function of the total stress 
intensity factor, K T . When presenting the results, the 
values for stress and stress intensity factors will be nor- 
malized by their values for a biaxial loading ratio value 
X = 0.5 such that: 

k t = (K t )/(K t ) 05 

(4) 

a = 

Thus, the normalized number of cycles for crack exten- 
sion, A, is rewritten as 


N = {\/K t ) A = (l/o) 4 (5) 

and can be used to estimate the effect of variations in 
stress results on durability. Similarly, the effect of varia- 
tions in stress results on the residual strength of the 
structure may be estimated by a normalized residual 
strength, 

RS = (1 /K t ) = (1/0). (6) 

To illustrate the effect of variations in stress results 
on durability and residual strength, A and RS are plotted 
versus {K T , a) in Fig. 4. As shown in Fig. 4, elevation 
of the stress results by 10%, 25%, and 50%, corresponds 
to reductions in A of 30%, 60%), and 80%, respectively. 



Kj y (5 

Figure 4. Normalized number of cycles to extend 
a crack. A, and normalized residual 
strength, RS , versus normalized stress 
intensity factor and normalized stress. 

Results and Discussion 

The results of a numerical study of the nonlinear re- 
sponse of a typical aircraft fuselage panel with and with- 
out longitudinal cracks, and subjected to combinations of 
internal pressure and axial compression loads are pre- 
sented. Results are presented first for a pristine structure, 
to demonstrate the effect of skin buckling on the fastener 
loads and the stress distribution in the fuselage panel, and 
to assess the impact of skin bucking on damage initia- 
tion. Then, results for fuselage panels with damage are 
presented to display the effect of skin bucking on the ten- 
dency for damage propagation. Results are presented for 
panels with damage along the skin-stringer attachment 
line in the form of failed fasteners or longitudinal skin 
cracks of various lengths. The results include contour 
plots of the circumferential stress and radial displace- 
ment in the fuselage skin and line plots of the stress in- 
tensity factors for cases with cracks. When presenting 
the stress intensity factor results, the components shown 
in Fig. 3 that are negligible for the current case are omit- 
ted from the line plots. 

Pristine Structure 

Nonlinear analyses were conducted for the fuselage 
panel with no damage. Solutions were obtained for an 
internal pressure load of 8 psi, and a range of axial load- 
ing values corresponding to biaxial loading ratio values 
X = 2.0 to -1.2. Contour plots of the fuselage-skin radial 
displacement for x = 1.5, 0.5, -0.5, and - 1 .0 are shown in 
Fig. 5. The case with x = 0.5 corresponds to the bulk- 
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| t — p-j— r 

0.0 Radial displacement, in. 0.15 


Frame Frame 



(b) X = 0.5 

(nominal pressure load) 



(c)x = -0.5 



(d) X - " 1 0 


Figure 5. Fuselage-skin radial displacement for 
biaxial loading ratio values of x = 1 -5, 
0.5, -0.5, and -1.0. 


head tension load for a nominal pressure load only. The 
radial displacement result for x = 0.5, shown in Fig. 5b. 
shows that the internal pressure on the skin deforms the 
skin radially outward, and the displacements are smaller 
where the skin is attached to the stiffening structure. The 
circumferential stiffness of the frames strongly resists 
the radial deformation, thus the skin’s radial displace- 
ment is the smallest near the frames. The stringers resist 
radial deformation because they are attached to the 
frames. The stringers provide some radial support to the 
skin through the bending stiffness of the stringers. The 
stringers bend along their length and deflect outward 
more than the frames, and the skin on each side of the 
stringer deflects outward more than the stringer. When 
the axial tension load is increased to a biaxial loading ra- 
tio value x = 1.5, the tension load tends to straighten the 
panel in the axial direction, and reduces the radial dis- 
placement between the frames, as shown in Fig. 5a. 
When the axial tension load is reversed to apply axial 
compression, solutions with negative biaxial loading ra- 
tio values are obtained. The skin radial displacement re- 
sults for the case with x = -0.5, shown in Fig. 5c. 
indicates that axial compression increases the outward 
deflection of the unsupported skin, and a displacement 
pattern with multiple axial half-waves between the 
frames starts to develop. The formation of additional 
half- waves is dominated by growth of the compres- 
sion-loaded bending boundary layer near each frame. 
When the axial compression load is increased beyond % 
= -0.6, the skin buckles into seven half waves between 
the frames, and one circumferential half wave between 
stringers. The skin radial displacement for x = -10, 
shown in Fig. 5d, displays a fully-developed buckled 
skin pattern that is symmetric with respect to each string- 
er and each frame. The symmetry in the response is at- 
tributed to the strong influence of the bending boundary 
layer on each side of the frames and the presence of the 
internal pressure load. Changes in the structural config- 
uration would likely influence the deformation shape and 
symmetry of the response. 

Panel cross sections ‘A’, k B' and ‘C’ are identified 
in Fig. 5. For each cross section, the deformed shape of 
the skin-stringer attachment area is shown amplified by 
a factor of 10 in Fig. 6. For the case with nominal pres- 
sure load, x = 0.5, the deformed shape of cross section 
4 A’, shown in Fig. 6a, indicates a small amount of out- 
ward deformation in the skin on each side of the stringer, 
and a small amount of twisting in the stringer deforma- 
tion due to the asymmetry of the stringer Z cross section. 
For the case with the postbuckled skin and x = -1.0, the 
deformed shape of cross section ‘B' shows the skin de- 
formed toward the stringer and bent over the stringer, 
with contact evident in the skin-stringer interface. At 
cross section ‘C’ of the postbuckled skin with x = -1.0, 
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30 


(a) x = 0.5, non-buckled skin, 
cross section ‘A 1 in Fig. 5b 



(b) x = - 1 -0, buckled skin, 
cross section 4 B' in Fig. 5d 



(c) x = - 1 0, buckled skin, 
cross section *C’ in Fig. 5d 


Figure 6. Deformed shape ( l Ox) of the stringer 
cross section for biaxial loading ratio 
values of X = 0.5 and - 1 .0. 

the deformed shape shown in Fig. 6c shows the skin 
pulled away from the stringer, which causes the asym- 
metric stringer to twist. The skin and stringer separate on 
one side of the fastener row (see Fig. 6c), and the bend- 
ing response of the skin is most severe in this region. 

Contour plots of the longitudinal stress on the inner 
and outer surface of the skin of the fuselage panel for val- 
ues of x = 0.5 and - 1 .0 are shown in Fig. 7. The area of 
the panel shown in the contour plots is one skin bay be- 
tween two frames and half a skin bay on each side of the 
center stringer. The results for the nominal pressure 
load, x = 0.5, indicate that the longitudinal stress in the 
skin is approximately 7 ksi in tension, and is nearly uni- 
form. For the case with % = - 1 .0, the buckled skin causes 
large gradients in the longitudinal stress. The buckled 
skin between the stringers has a middle surface stress of 
approximately -5 ksi, but the bending deformation caus- 
es surface stresses ranging from -23 ksi to 13 ksi. The 
buckled skin has a lower longitudinal stiffness than the 
unbuckled skin. Thus, a higher percentage of the axial 
load is supported by the stringers and the skin near the 
stringers. For the case with x = - TO, the skin near the 
stringers has a longitudinal stress of approximately 
-14 ksi. 


-30 Longitudinal stress, ksi 


Outer surface 



(a) x ~ 0.5 (nominal pressure load) 
unbuckled skin 


Outer surface 



(b) x = - 1 0, buckled skin 


Figure 7. Fuselage-skin longitudinal stress on 
inner and outer surfaces for biaxial 
loading ratio values of x = 0.5 and - 1 .0. 
(skin region between frames and over 
center stringer). 


7 

American Institute of Aeronautics and Astronautics 






The stringers are subjected to axial loads that are 
applied to the ends on the panel, and to bending loads. 
The bending loads result from the internal pressure on 
the panel skin, and the radial constraint imposed by the 
frames. The stresses in the stringers are predominately 
longitudinal, and the maximum stresses are located in the 
cap of the stringer, i.e., away from the skin-to-stringer 
connection. For the nominal pressure load, x - 0.5, the 
longitudinal stress in the cap of the stringer is approxi- 
mately 8 ksi at the frame location, and -1 ksi midway be- 
tween frames. The longitudinal stress in the attachment 
flange of the stringer is approximately 7 ksi. For the 
compression loaded case with % = -1 0, the longitudinal 
stress in the cap of the stringer is approximately -8 ksi at 
the frame location, and -29 ksi midway between the 
frames. The longitudinal stress in the attachment flange 
of the stringer is approximately - 1 5 ksi. Thus, the bend- 
ing load on the stringer is greater when the panel is load- 
ed in compression and the skin is buckled. 

Contour plots of the circumferential stress on the 
inner and outer surface of the skin of the fuselage panel 
for values of % = 0.5 and -1.0 are shown in Fig. 8. The 
area of the panel shown in the contour plots is one skin 
bay between two frames and half a skin bay on each side 
of the center stringer. The results lor the nominal pres- 
sure load, x ~ 0-5, indicate that the circumferential stress- 
es in the skin are nearly uniform, and are smaller closer 
to the frames. Along the line of skin-stringer attachment, 
the skin bending shown in Fig. 6a causes an elevation in 
the circumferential stress on the inner surface and a re- 
duction in the circumferential stress on the outer surface. 
For the case with x = - 1 -0, the buckled skin causes large 
gradients in the circumferential stress. In regions where 
the skin buckles outward, the circumferential stresses are 
elevated. In regions where the skin buckles inward, the 
circumferential stresses are reduced. The largest stresses 
are located along the line of skin-stringer attachment. At 
points ‘D' and *E’ on the panel skin, identified in Fig. 8b, 
there are significant bending stresses associated with the 
skin bending shown previously in Figs. 6b and 6c. The 
variation in circumferential stress at skin locations ‘D’ 
and 4 E’ for biaxial loading ratio values of x = 1 .5 to - 1 .2 
is shown in Fig. 9. In this figure, inner-, middle-, and 
outer-surface stresses are presented, and the results are 
normalized by the middle-surface stress for x = 0.5. For 
values of x = 15 to -0.5, the middle-surface circumfer- 
ential stresses are insensitive to the changes in the axial 
load, and the inner- and outer-surface stresses indicate 
that axial compression causes a linear increase in the 
bending stress that increases the inner-surface tensile 
stress, and decreases the outer-surface tensile stress. For 
values of X = -0.5 to - 1 .2, skin buckling causes nonlinear 
changes in the middle-surface and bending stresses. The 
tensile circumferential stresses on the outer surface at lo- 


-5 Circumferential stress, ksi 30 

Outer surface 


x Inner surface 


(a) x = 0.5 (nominal pressure load) 
unbuckled skin 

Outer surface 


x Inner surface 


(b) x = - 1 .0, buckled skin 

Figure 8. Fuselage-skin circumferential stress on 
inner and outer surfaces for biaxial 
loading ratio values of x = 0.5 and - 1 .0. 
(skin region between frames and over 
center stringer). 
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Figure 9. Circumferential stress at skin locations 
‘D’ and ‘E’ identified in Fig. 8b. for 
biaxial loading ratio values of % = 1.5 
to - 1 .2; inner-, middle-, and 
outer-surface stresses normalized by 
the middle- surface stress for x = 0.5. 

cation ‘D\ and on the inner surface at location *E\ are 
more than two times the values for the pressure-only 
case. In addition, the middle- surface tensile stress at lo- 
cation ‘E’ increases by approximately 30 percent. The 
large stress values located in the skin-stringer attachment 
region will increase the likelihood of damage initiation 
and propagation in this region. Large surface stresses 
will promote initiation of a surface crack, while a large 
middle- surface stress will promote growth of a 
through-crack. 

The maximum forces in the fasteners that connect 
the skin to the center stringer are shown in Table 1 for bi- 
axial loading ratio values x = 0.5 and -1.0. The axial 
shear forces are oriented along the stringer longitudinal 
axis, while the side shear forces are perpendicular to the 
stringer longitudinal axis. The bending moment that is 
reported is the moment about the stringer longitudinal 
axis. The magnitude of the bending moments are small 
compared to the forces, but, considering that the moment 
arm of the fastener is also small, the moments may be 
significant. When evaluating the maximum forces, the 
fasteners in the center two-thirds of the stringer length 
between the frames were considered. The results in 
Table 1 indicate that the maximum fastener loads are 
along the stringer longitudinal axis. This condition will 
be the case, in general, since the stringer is much stiffer 
along its longitudinal axis than in any other direction. 
Compared to the pressure-only case, x = 0*5, cases with 


Table 1. Maximum fastener forces along the center 


stringer in a pristine panel 


Biaxial 
loading 
ratio, x 

Tension, 

lb. 

Axial 

shear, 

lb. 

Side 

shear, 

lb. 

Bending 

moment, 

in-lb. 

0.5 

8 

40 

5 

2 

-1.0 

50 

1 14 

58 

8 


pressure plus axial tension, x > 0.5, have smaller fastener 
forces. When axial compression is applied and the skin 
buckles, as is the case for x = -L0, all of the fastener 
loads become significantly larger. The maximum tensile 
load in the fasteners occurs at locations where the skin 
deforms radially outward on each side of the stringer. 

If the panel skin has buckled, the distribution of fas- 
tener forces is quite complex. The short- wavelength 
buckling deformation has only three fasteners per axial 
half-wave. The large displacement gradients result in 
fastener forces that are rather erratic. Also, the maxi- 
mum fastener force for each component may refer to sev- 
eral different fasteners, rather than a single fastener. The 
erratic nature and complexity of the response may make 
it difficult to determine a critical location, and make it 
difficult to define a refined local model to examine the 
fastener detail more closely. 

Three Failed Fasteners 

One type of damage that may occur in a fuselage 
panel is failure of the skin-to-stringer attachment. In the 
results described above, it was noted that buckling of the 
panel skin increases the tensile forces at the skin-stringer 
attachment at locations where the skin deforms radially 
outward on each side of the stringer. These higher at- 
tachment loads increase the likelihood of fastener failure 
by fatigue, or through the formation of small cracks at 
the fastener holes which would promote pull-through of 
the fasteners. To examine the effect of fastener failure on 
the panel’s response, three consecutive fasteners in prox- 
imity to point ‘E’ in Fig. 8b were removed from the finite 
clement model, and nonlinear analyses were conducted. 
Contour plots of the radial displacement of the fuselage 
skin with three fasteners removed are shown in Fig. 10 
for biaxial loading ratio values of x = 0.5 and - 1 .0. For 
the nominal pressure load, x = 0.5, the displacement is 
only slightly different from the displacement of the pris- 
tine panel shown in Fig. 5b. For the case with axial com- 
pressive load, x = - LO, the panel skin is buckled, and the 
maximum radial displacement occurs in the skin where 
the fasteners were removed, or location ‘F* indicated in 
Fig. 10b. The deformed shape has a single half wave in 
the circumferential direction that extends across two skin 
bays, and the magnitude of the displacement near the re- 
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0.0 


Radial displacement, in. 


0.15 


Frame Frame 



(a) x = 0-5 

(nominal pressure load) 



(b)x = -1.0 


Figure 10. Fuselage-skin with three fasteners 
removed: radial displacement for 
biaxial loading ratio values of X = 0.5 
and -1.0. 


moved fasteners is 20% larger than the maximum dis- 
placement of the pristine panel shown in Fig. 5d. 

The variation in circumferential stress at skin loca- 
tion 4 F’ for biaxial loading ratio values of x = 1 .5 to - 1 .2 
is shown in Fig. 1 1. Inner-, middle-, and outer-surface 
stresses are presented in this figure, and the results are 
normalized by the middle-surface stress for the case with 
no fasteners removed, and x = 0.5. For values of x = 1 .5 
to 0.0, the normalized stresses are approximately equal to 
1 .0. The stresses are insensitive to the changes in the ax- 
ial load, and there is very little bending in the skin. When 
the axial load becomes compressive, x = 0.0 to -1 .2, the 
skin buckles and the tensile circumferential stress in- 
creases dramatically, more than doubling in magnitude at 
the middle and outer surfaces. The large stress values lo- 
cated in the skin region near the removed fasteners will 
increase the likelihood of the formation and propagation 
of skin cracks in this region. Refer to Fig. 4 to estimate 
the effect of variations in stress results on the durability 
and residual strength of the panel. 



Figure 1 1. Circumferential stress at skin location 
4 F’ identified in Fig. 10b. for biaxial 
loading ratio values of x = 1.5 to - 1 .2; 
inner-, middle-, and outer-surface 
stresses normalized by the 
middle-surface stress for the case with 
no fasteners removed, and x = 0.5. 

Centered 1 -In. -Long Longitudinal Crack 

The finite element model was modified to include a 
1 -in. -long longitudinal crack in the panel skin. The crack 
was located midway between frames, centered on x=0 
in Fig. 1 , and adjacent to the line of skin-stringer attach- 
ment. Referring to the stringer cross section shown in 
Fig. 6, the crack was located 0.25 in. to the left of the fas- 
tener location. Nonlinear analyses were conducted for 
an internal pressure load of 8 psi, and a range of axial 
load values corresponding to biaxial loading ratio values 
X = 2.0 to - 1 .2. The 1 -in. -long crack did not influence the 
panel’s overall response. When the panel was loaded in 
compression, the skin buckled into the same deformed 
shape as observed for the pristine panel. The fuse- 
lage-skin radial displacement in the region near the 
crack, and with x = -1.0, is shown in Fig. 12. The dis- 
placement results indicate that the skin nearest to the 
crack buckles radially inward, which, for the pristine 
panel, was shown in Fig. 9 to reduce the middle-surface 
circumferential stress locally. Radial displacement of 
buckled skin with 1 -in. -long longitudinal crack midway 
between frames (x = 0.0) and adjacent to the center 
stringer. 

Stress intensity factors for the 1 -in. -long longitudi- 
nal crack, and biaxial loading ratio values of x = 2 to -1 .2 
are shown in Fig. 1 3. The stress intensity factors are nor- 
malized by K t = 15 ksivin. for x = 0.5. Results are 
shown for Kj , the symmetric membrane component, Kj , 
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x = 0 


Figure 12. Radial displacement of buckled skin 
with 1 -in. -long longitudinal crack 
midway between frames (x = 0) and 
adjacent to the center stringer, x = - 1 .0. 



K 


Figure 13. Stress intensity factors fora 1-in. -long 
longitudinal crack located midway 
between frames (x = 0) and adjacent to 
the center stringer, for biaxial loading 
ratio values of Y = 2 to -1 .2; normalized 
by Kj — 15 ksiVin. for x = 0.5. 

and the symmetric bending component, k h For biaxial 
loading ratio values X = 2.0 to -0.5, the crack-tip re- 
sponse is dominated by Kj , and the response is not very 
sensitive to variations in the biaxial loading ratio. When 
an axial compression load is applied and the skin buckles 
(X < -0.6), the stress intensity lactor Ay is reduced in mag- 
nitude, and the bending component is increased in mag- 


nitude. When the skin buckles, the crack-tip response at 
opposite ends of the crack are similar, but not identical, 
because the center of the crack and the skin deformation 
pattern are not aligned in the axial direction. The buck- 
led skin causes high gradients in the response, and small 
changes in the location of the crack-tip relative to the 
buckling deformation can influence the crack-tip re- 
sponse. The I -in. -long crack has no significant affect on 
the fastener forces or stringer stresses compared to the 
results for the pristine panel. 

Off-Center l-In.-Long Longitudinal Crack 

In the previous case, the 1 -in. -long crack was locat- 
ed midway between frames, centered on x = 0. In the 
current case, the 1 -in. -long crack is still adjacent to the 
line of skin-stringer attachment, but is shifted axially, 
and centered on x = -2.5 in Fig. 1. Again, the I. in. -long 
crack does not influence the panel's overall response. 
The fuselage-skin radial displacement in the region near 
the crack, and with x = -1.0 is shown in Fig. 14. The dis- 
placement results show that the skin nearest to the crack 
buckles radially outward, which, for the pristine panel, 
was shown in Fig. 9 to increase the middle-surface cir- 
cumferential stress locally. 

Stress intensity factors for the 1 -in. -long longitudi- 
nal crack centered on x = -2.5. and biaxial loading ratio 
values ol x = 2 to -1.2, are shown in Fig. 1 5. The stress 
intensity factors are normalized by K r - 15 ksiVin. for 
X — 0.5. Results are shown for Kj , the symmetric mem- 
brane component, Kj , and the asymmetric bending com- 
ponent, k 2 . For biaxial loading ratio values x = 2.0 to 
-0.5, the crack-tip response is dominated by K h and the 
response is not very sensitive to variations in the biaxial 
loading ratio. When axial compression load is applied 
and thej>kin buckles (x < -0.6), the stress intensity factors 
Kj and k 2 increase in magnitude. The symmetric mem- 
brane component, K h for x = - 1 .2 is 40 % larger than Kj 
for the pressure-only case, x = 0.5. The asymmetric 
bending component, f 2 , which is insignificant for 
X = 0.5, increases in magnitude by over 1.200% for 
X = -1.2, and is 80% larger than Ay for the pressure-only 
case, x = 0.5. Thus, the magnitude and mode-mixity of 
the crack-tip stress intensity factors are strongly influ- 
ence by skin buckling. By comparing these results with 
the results with the crack centered on x = 0. it is evident 
that the effect that skin buckling has on the crack-tip re- 
sponse of a short crack is dependent on the crack location 
relative to the buckling deformation. Refer to Fig. 4 to 
estimate the effect of variations in the total stress intensi- 
ty factor, A 7 ; on the durability and residual strength of the 
panel. 
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0.0 

Radial displacement, in. 

0.15 


i -in. -long crack 


Fasteners for 
center stringer 


x = -2.5 

Figure 14. Radial displacement of buckled skin 
with a 1 -in. -long longitudinal crack 
centered on x = -2.5 in. and adjacent to 
the center stringer, X = ' 1 0. 



Figure 15. Stress intensity factors for a I -in. -long 
longitudinal crack centered on x = -2.5 
in. and adjacent to the center stringer, 
for biaxial loading ratio values of x f 2 
to - 1 .2; normalized by K r = 1 5 ksiVin. 
for % = 0.5. 

Centered 4-In.-Long Longitudinal Crack 

The finite element model was modified to include a 
4-in. -long longitudinal crack in the panel skin. The crack 
was located midway between frames, centered on x = 0 
in Fig. 1 , and adjacent to the line of skin-stringer attach- 
ment. Nonlinear analyses were conducted for an internal 
pressure load of 8 psi, and a range of axial load values 


corresponding to biaxial loading ratio values x = 2.0 to 
-1.2. The 4-in. -long crack is large enough to influence 
the panel’s overall response. Contour plots of the fuse- 
lage-skin radial displacement for biaxial loading ratio 
values x = 0.5 and -1.0 are shown in Fig. 16. The dis- 
placement results indicate that the presence ot the 
4-in. -long crack influences the radial displacement in 
one skin bay on each side of the center stringer. The in- 
fluence is not confined to the skin bay on the side of the 
stringer where the crack is located, because the crack un- 
loads the circumferential tension load in the skin, and the 
stringer is not stiff enough in the circumferential direc- 
tion to prevent the adjacent skin bay from also unloading. 

A contour plot of the middle-surface longitudinal 
stress in the fuselage-skin is shown in Fig. 17 for the 
nominal pressure load, x = 0.5. The results shown in 
Fig. 17 suggest that the influence of the 4-in. -long crack 
is contained within the two-skin-bay area in the center of 
the panel. The results in Fig. 17 also indicate that there 
is a small region of the skin on each side of the crack that 
has a compressive longitudinal stress. These compres- 
sive stresses are like the compressive stresses that devel- 
op in regions near the edge of the crack in a 
center-cracked flat plate subjected to a tension load. The 
stresses in the stringer are influenced by the crack in the 
skin in the following manner. The skin on one side of the 
4- in-long crack is attached to the stringer. Along the 
length of the crack, the compressive longitudinal stress 
in the skin applies a compressive load to the stringer at- 
tachment flange. In addition, the remote circumferential 
stress in the skin causes the crack to open, which imparls 
a circumferential side load on the stringer. 

For cases with all values ot the biaxial loading ra- 
tio, the radial displacement is larger in the two skin bays 
adjacent to the crack, and the shape of the buckling de- 
formation is different from the deformation in the re- 
mainder of the panel. The internal pressure causes 
outward bulging of the skin near the crack, and these 
bulging deformations are magnified when the panel is 
subjected to compressive loads. The deformations asso- 
ciated with the 4-in. -long crack dominate the local panel 
response and skin buckling deformation. This behavior 
was not the case for shorter cracks, where the skin buck- 
ling deformation dominated the panel response, and the 
crack behaved more like a local detail or inclusion. The 
deformed shape (3x magnification) of the center stringer 
near the 4-in.-long crack for biaxial loading ratio values 
of x = 0.5 and - 1 .0 are shown in Fig. 1 8. The deformed 
shapes have significant displacements in the skin, but 
distortion of the stiffener cross-section appears to be 
minimal. The results of the analyses indicate that the 
stringer did not yield or collapse, and was able to support 
the additional loads developed by the crack. 
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Figure 17. Middle-surface longitudinal stress in 
the fuselage-skin with a 4-in. -long 
longitudinal crack located midway 
between frames (x = 0.0) and adjacent 
to the center stringer, for biaxial 
loading ratio value of % = 0.5. 


Radial displacement of fuselage-skin 
with a 4-in. -long longitudinal crack 
located midway between frames 
(x = 0.0) and adjacent to the center 
stringer, for biaxial loading ratio values 
of X = 0.5 and -1.0 


Figure 18. Deformed shape (3x) of the center 
stringer near a 4-in. -long crack for 
biaxial loading ratio values of y = < 
and -1.0. 
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The maximum forces in the fasteners that connect 
the skin to the center stringer are reported in Table 2 for 
biaxial loading ratio values x = 0.5 and -1.0. The maxi- 
mum fastener loads are considerably larger than the val- 
ues for the pristine panel reported in Table 1 . Compared 
to the pressure-only case, X = 0.5, cases with pressure 
plus axial tension, x > 0-5, have smaller fastener forces. 
When axial compression is applied and the skin buckles, 
as is the case for x = -1.0, all of the fastener loads be- 
come significantly larger. 


Table 2. Maximum fastener forces along the center 
stringer in a panel with a 4-in. -long 
longitudinal crack. 


Biaxial 
loading 
ratio, X 

Tension, 

lb. 

Axial 

shear, 

lb. 

Side 

shear, 

lb. 

Bending 

moment, 

in-lb. 

0.5 

18 

58 

55 

3 

-1.0 

81 

328 

115 

8 


Stress intensity factors for the 4-in. -long longitudi- 
nal crack are shown in Fig. 19 for biaxial loading ratio 
values of x = 2 to -1.2. The stress intensity factors are 
normalized by Kj = 36 ksiVin. lor x = 0.5. Results are 
shown for K T < the symmetric and asymmetric membrane 
components, Kj and Kjj , respectively, and the asymmet- 
ric bending component, k 2 . For the pressure-only case, 
X — 0.5, and cases with pressure plus axial tension. 



Figure 19. Stress intensity factors for a 4-in. -long 
longitudinal crack centered between 
frames (x = 0) and adjacent to the 
center stringer, for biaxial loading ratio 
values of x = 2 to - 1 .2; normalized by 
K r = 36 ksh/in. for x = 0.5. 


X > 0.5, the crack-tip response is dominated by K h and 
the response is not very sensitive to variations in the bi- 
axial loading ratio. When axial compression load is ap- 
plied, X < 0.5, the stress intensity factors increase in 
magnitude in a manner that is typical of a limit-load re- 
sponse, rather than a bifurcation buckling response. That 
is, the bulging deformation near the crack develops grad- 
ually with increasing compression load, rather than 
changing suddenly when the skin buckles. For the max- 
imum compression load considered, x = -1*2, Kj is 70% 
larger than Kj for the pressure-only case. Similarly, k 2 
and Mu are 270% and 460% larger, respectively, for 
X = - 1.2, than their respective values for a pressure load 
only, x = 0.5. Refer to Fig. 4 to estimate the^eflect ot 
variations in the total stress intensity factor. K r . on the 
durability and residual strength of the panel. 

Centered lO-In.-Long Longitudinal Crack 

The final crack configuration considered in the 
present study is a lO-in.-long longitudinal crack located 
midway between frames, centered on x = 0 in Fig. 1 , and 
adjacent to the line of skin-stringer attachment. Nonlin- 
ear analyses were conducted and the panel's response 
was quite similar to that obtained for the 4-in. long crack. 
Contour plots of the fuselage-skin radial displacement 
for biaxial loading ratio values x = 0-5 and -1.0 are 
shown in Fig. 20. The displacement results indicate that 
the presence of the longer crack influences the radial dis- 
placement over a larger portion of the panel. 

A contour plot of the middle-surface longitudinal 
stress in the fuselage-skin is shown in Fig. 21 for the 
nominal pressure load, x = 0*5. Middle-surface stress in 
fuselage-skin with lO-in.-long longitudinal crack mid- 
way between frames (x = 0.0) and adjacent to the center 
stringer. The results shown in Fig. 21 suggest that the in 
fluence of the lO-in.-long crack extends over the entire 
length and width of the panel. The results in Fig. 21 also 
indicate that there is a region of the skin that has a com- 
pressive longitudinal stress, on the side of the crack that 
is attached to the stringer. Recall that the results for the 
case with a 4-in. -long crack, shown in Fig. 17, indicated 
compressive longitudinal stresses in the skin on both 
sides of the crack. When the crack is longer, the internal 
pressure load on the skin causes larger outward deforma- 
tion on the unsupported side of the crack, and this defor- 
mation induces a nonlinear tensile stress along the 
unsupported crack edge. This tensile longitudinal stress 
along the unsupported crack edge acts to increase the 
compressive longitudinal stress on the side of the crack 
that is supported by the stringer. The stresses in the 
stringer are influenced by the crack in the skin in the fol- 
lowing manner. The skin on one side of the 10-in-long 
crack is attached to the stringer. Along the length of the 
crack, the compressive longitudinal stress in the skin that 
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Figure 20. Radial displacement of a fuselage skin 
with a 10-in. -long longitudinal crack 
located midway between frames 
(x = 0.0) and adjacent to the center 
stringer, for biaxial loading ratio values 
of x — 0.5 and -1.0. 



Figure 21. Middle-surface longitudinal stress in a 
fuselage skin with a 10-in. -long 
longitudinal crack located midway 
between frames (x = 0.0) and adjacent 
to the center stringer, for biaxial 
loading ratio value of X = 0.5. 

is supported by the stringer applies a compressive load to 
the stringer attachment flange. In addition, the remote 
circumferential stress in the skin causes the crack to 
open, which imparts a circumferential side load on the 
stringer. 

The deformed shape (3x magnification) of the cen- 
ter stringer near the 10-in. -long crack is shown in Fig. 22 
for biaxial loading ratio values of x = 0.5 and - 1 .0. The 
deformed shapes indicate that there are significant dis- 
placements in the skin, and the stringer cross-section ap- 
pears to be rolling for the biaxial loading ratio value of 
X = -10. A review of the analysis results indicates that 
the center stringer exhibits a limit load behavior at a bi- 
axial loading ratio value of x = - 1 0, and the stress levels 
in several locations of the model exceed the yield stress 
for the 2024-T3 aluminum alloy. 

The maximum forces in the fasteners that connect 
the skin to the center stringer are shown in Table 3 for bi- 
axial loading ratio values x = 0.5 and -1.0. The maxi- 
mum fastener loads are considerably larger than the 
values for the panel with the 4-in. long crack reported in 
Table 2. Compared to the pressure-only case, x = 0.5, 
cases with pressure plus axial tension, x > 0.5, have 
smaller fastener forces. When axial compression is ap- 
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Figure 22. Deformed shape (3x) of the center 
stringer near a 10-in. -long crack for 
biaxial loading ratio values of X = 0*5 
and - 1 .0. 


Table 3. Maximum fastener forces along the center 
stringer in a panel with a 10-in. -long 
longitudinal crack 


Biaxial 
loading 
ratio, X 

rp Axial 

Te ™ on ' shear, 

lb - lb. 

Side 

shear, 

lb. 

Bending 

moment, 

in-lb. 

0.5 

64 

237 

109 

7 

-1.0 

139 

641 

283 

12 


plied and the skin buckles, as is the case for X - -1-0, all 
of the fastener loads become significantly larger. 

Stress intensity factors for the 10-in. -long longitu- 
dinal crack, and biaxial loading ratio values of % = 2 to 
*1 .2 are shown in Fig. 23. The stress intensity factors are 
normalized by K T =12 ksiVin. for x = 0.5. Results are 
shown for K T> the symmetric and asymmetric membrane 
components, ~Kf and K fI , respectively, and the asymmet- 
ric bending component, I?* The crack-tip response is 
dominated by Kj for all values of load, and for cases with 
pressure plus axial tension, x > 0.5, K } is not very sensi- 
tive to variations in the biaxial loading ratio. When an 
axial compression load is applied, x < 0.5, the stress in- 
tensity factors increase in magnitude, although the per- 
centage increases are less than for the case with the 
4-in.-long crack. For the maximum compression load 
considered, % = - 1 .2, K I is 56% larger than Kj for the 
pressure-only case. Similarly, £? an ^ K/i are 58% and 
120% larger, respectively, for x = - 1 -2, than their respec- 
tive values for pressure-only load, X = 0.5. Refer to 



Figure 23. Stress intensity factors for a 1 0-in.-long 
longitudinal crack centered between 
frames (x = 0) and adjacent to the 
center stringer, for biaxial loading ratio 
values of x = 2 to - 1 .2; normalized by 
K r =12 ksiVin. for x = 0.5. 

Fig. 4 to estimate the effect of variations in the total 
stress intensity factor, K T , on the durability and residual 
strength of the panel. 

Concluding Remarks 

The results of a numerical study to assess the effect 
of skin buckling on the internal load distribution in a 
stiffened fuselage panel, with and without a longitudinal 
skin crack, are presented. In addition, the impact ol 
changes in the internal loads on the fatigue life and resid- 
ual strength of a fuselage panel is assessed. The results 
indicate that nonlinear analyses of the stiffened-shell 
model can provide accurate predictions of the geomet- 
ric-nonlinear response of the buckled skin, cross section 
deformation of the stiffening components, and 
skin-stringer attachment forces associated with discrete 
fasteners. The numerical results indicate that compres- 
sion loads and skin buckling can have a significant effect 
on the circumferential stress in the skin, and fastener 
loads, which will influence damage initiation. Compres- 
sion loads and skin buckling have a comparable effect on 
stress intensity factors for cases with cracks, which will 
influence damage propagation rates and the residual 
strength of the panel. 

The results from analyses of a pristine panel indi- 
cate that skin buckling creates large gradients in the radi- 
al displacement, longitudinal stresses, and 
circumferential stresses in the fuselage-panel skin. The 
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largest tensile stresses in the panel are located where the 
skin buckles radially outward on each side of the string- 
er. The fastener tension loads, normal to the skin surface, 
are also the largest in these areas. The increase in the cir- 
cumferential stress near the stringer due to skin buckling 
is quantified, and the associated reduction in fatigue life 
is estimated to be significant. Results for analyses with 
fasteners removed indicated that fastener failure, by fa- 
tigue or pull through, would have little effect on the re- 
sponse of a panel under axial tensile load. But if a 
compression load causes the skin to buckle, then re- 
moved fasteners can increase the circumferential stress 
in the skin, and reduce the fatigue life. 

The numerical results for cases with longitudinal 
cracks indicate that compressive loads and skin buckling 
generally increases the crack-lip stress intensity factors, 
and thus, reduces the fatigue life and residual strength of 
the panel. For cases with short cracks, the global panel 
response was unaffected by the crack, and the crack-tip 
response was sensitive to the crack location relative to 
the buckled shape of the panel. For cases with longer 
cracks, axial compression loads couple with the radial 
deformations near the crack to create a limit-load type of 
response. The region of the panel that is influenced by 
the crack is larger. Also, axial compressive loads and 
side loads on the stringer adjacent to the crack are in- 
creased by the presence of a crack, and stringer strength 
or stability may become an issue. 
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